Computer simulations of colloidal transport on a patterned magnetic substrate 
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We study the transport of paramagnetic colloidal particles on a patterned magnetic substrate 
with kinetic Monte Carlo and Brownian dynamics computer simulations. The planar substrate is 
decorated with point dipoles in either parallel or zigzag stripe arrangements and exposed to an 
additional external magnetic field that oscillates in time. For the case of parallel stripes we find 
that the magnitude and direction of the particle current is controlled by the tilt angle of the external 
magnetic field. The effect is reliably obtained in a wide range of ratios between temperature and 
magnetic permeability. Particle transport is achieved only when the period of oscillation of the 
external field is greater than a critical value. For the case of zigzag stripes a current is obtained 
using an oscillating external field normal to the substrate. In this case, transport is only possible in 
the vertex of the zigzag, giving rise to a narrow stream of particles. The magnitude and direction of 
the particle current are found to be controlled by a combination of the zigzag angle and the distance 
of the colloids from the substrate. Metropolis Monte Carlo and Brownian dynamics simulations 
predict results that are in good agreement with each other. Using kinetic Monte Carlo we find that 
at high density the particle transport is hindered by jamming. 
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I. INTRODUCTION 

Manipulation and transport of magnetic particles at 
the nanometer and micrometer scales are important tech- 
nological processes for biological [1] and biomedical ap- 
phcations [2]. Magnetic microspheres are routinely used 
as markers for cells or larger molecules. Micron-sized 
magnetic beads were attached to the free end of viral 
DNA 3 for the direct visualization in a bright field mi- 
croscope of a DNA-packaging process. Molday et al. Q 
achieved separation of red blood cells and lymphoid cells 
using magnetic particles chemically bound to antibodies. 
Pamme and Wilhelm [5| achieved continuos cell sorting of 
magnetically labeled cells via free-flow magnetophoresis. 
Magnetic particles can also be tagged with a fluorescent 
dye and followed in real time using microscopy [6|, and 
their surface can be functionalized to selectively bind to 
specific targets in a solution, allowing for selective sepa- 
ration using magnetic fields 0- 

Since a magnetic particle in a homogenous magnetic 
field cannot have a net translational motion, a variety of 
methods have been developed for the generation of inho- 
mogeneous magnetic fields that can be used to manipu- 
late and transport magnetic particles. Deng et al. [H and 
Lee et al. [9] used lithography to fabricate circuits that 
carry electrical currents in order to generate the inhomo- 
geneous magnetic field that is necessary for manipulating 
the magnetic particles. Another approach relies on the 
deposition of discrete ferromagnetic elements. In particu- 
lar, Yellen et al. [lo| printed cylindrical magnetic islands 
on a substrate to drive the assembly of colloidal particles. 
The application of a rotating external magnetic field al- 
lows to control and transport small super-paramagnetic 
particles [l2|,|T3|. Gunnarsson et al. [ll| deposited ellipti- 
cal magnetic elements and showed the ability to transport 
paramagnetic particles along the ellipses by the applica- 
tion of a rotating magnetic field. Another technique uses 



ferrite garnet films [l4] , which show patterns of alternat- 
ing magnetization, and an external field that oscillates in 
time in order to control the transport of colloidal parti- 
cles 

[ii-Iilj- The pattern on a garnet film forms spon- 
taneously and is not easily controllable. On the other 
hand, the deposition of small magnetic islands [lo|, HH, 
opens up the possibility of creating structured magnetic 
substrates with full control of the deposition pattern. 

In this hght, we carry out both kinetic Metropolis 
Monte Carlo and Brownian Dynamics computer simu- 
lations to study the behavior of paramagnetic colloidal 
particles on a substrate of discrete magnetic dipoles ar- 
ranged in two specific patterns, namely parallel and 
zigzag stripes. The patterns were inspired by those 
found in garnet films [14-16], but could, in principle, 
also be fabricated by one of the above deposition tech- 
niques [io|,[Iii. The field generated by a garnet film [l5| 
is quantitatively different from the one produced by an 
array of discrete dipoles. Nevertheless, we show that the 
differences are small and that our model shows a par- 
ticle transport behavior similar to the one found exper- 
imentally on garnet films. In particular, we find that 
the random Brownian motion of the colloidal particles 
is turned into a deterministic motion by an oscillating 
external magnetic field. We analyze the conditions that 
enable the deterministic motion and hence the controlled 
transport of the paramagnetic colloidal particles. 

Brownian dynamics is based on the equations of mo- 
tion for overdamped particles without hydrodynamic in- 
teractions, on the other hand Monte Carlo reproduces 
the correct dynamics only under certain conditions [itI - 
[22| . The advantage of Monte Carlo is its higher compu- 
tational efficiency with respect to Brownian dynamics. 
Therefore, we explicitly compare the results of Monte 
Carlo and Brownian dynamics at low density, and carry 
out only Monte Carlo simulations in the computationally 
demanding high density regime. 
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Figure 1: (Color online) Sketch of the model. The colloidal 
particles (spheres) are suspended at a distance Zcoii over a 
pattern of discrete point dipoles (arrows). The dipoles are 
oriented along the x-axis with alternating magnetization m = 
(zbmo, 0, 0). The wavelength of the repeating pattern is A. 
The particles in the figures are shown with diameter 2cr to 
help visualization, (a) Parallel stripes pattern, (b) Zigzag 
pattern with the zigzag angle Oz as a control parameter. 



The paper is organized as follows: In Sec. [Ill we summa- 
rize the model used and describe the simulation details. 
In particular, in Sec. Ill Al we define the particle-particle 
and particle-substrate interactions, in Sec. IIIBI we dis- 
cuss the energy landscape of the model, and in Sec. Ill CI 
we give the simulation details. In Sec. HIH we show and 
discuss the results for both parallel and zigzag stripes. In 
SecIVlwe give some concluding remarks. In appendix [Al 
we discuss similarities and differences between the field 
produced by our model and that produced by a Garnet 
film. 



II. MODEL AND METHOD 



A. Definition of the interactions 



dependent and spatially homogeneous external magnetic 
field Hext(^) = H^^^ sin(27rt/ro) is applied to the sys- 
tem. Here H^^"" = {H^,Hy,H^) is the amplitude of the 
external field, tq is its oscillation period, and t is the 
time. 

A colloidal fluid of N paramagnetic spheres with hard- 
core diameter a lies suspended at a distance Zcoii from 
the patterned substrate and is constrained to move in 
the x — y plane only. The total magnetic field exerted on 
a paramagnetic particle i at position = (x^, yi^ Zcow) is 
the sum of the external field and the substrate field 



H(r„t) = Hsub(r.) + HextW 



(2) 



Hence, a dipole moment = xH(ri,t), is induced in 
the paramagnetic particle i with susceptibility x- 

The interaction energy between the dipoles in the sub- 
strate is constant in time, therefore the relevant energy 
of our model is the sum of three contributions: first the 
hard-core interaction between the particles, second the 
interaction between the particles' (induced) dipole mo- 
ments and the total magnetic field H(x, Zcoii, 
and third the dipole-dipole interaction between the par- 
ticles. The total energy can therefore be written as 

/3Utot{t) = ^$(n,) - /3/i.x^H(r.,t)2 

i<j i 



(3) 



^ [3 H(r„ t) . e,, H(r,, • e,, - H(r„ t) • H(r„ ^)] , 



i<j 



where the hard-sphere potential ^{ri^ 



oo if ra < a and 



Ir^ — r. l the distance between col- 



loidal particles i and j, e^^ 



,yi, Zcou) and 



— i^j^Vj^ ^^coii)- The field H is defined by equation (|2]), /is 
is the magnetic permeability of the solvent and /3 = l/ksT^ 
where /cs is the Boltzmann constant and T is the tempera- 
ture. 



We study a fixed lattice of size N^x Ny^ of point mag- 
netic dipoles min lying in the x-y plane with components 
(zbmo, 0, 0). The lattice sites are enumerated by the pair 
of integers (/,n), where / refers to the x-direction and n 
to the ^/-direction. The x-component is +mo for dipoles 
sitting at an odd / position and —mo for dipoles sitting at 
an even / position. The dipole moments have all the same 
magnitude mo and form a pattern of parallel (Fig. [TJi) or 
zigzag stripes (Fig.[T]3). The wavelength of the repeating 
pattern in the x-direction is denoted by A. The sepa- 
ration distance between point dipoles in the ^/-direction 
is /S.y and the separation distance between point dipoles 
in the x-direction is Ax = A/2. The zigzag pattern is 
characterized by the angle Oz^ as shown in Fig. [Ud. 

The substrate generates a magnetic field 

where Yin is the distance between the dipole (/,n), and 
the space point r = (x^y^z). In addition, a time- 



B. Analysis of the energy landscape 

As shown by Eq. (j3j), the particle-particle interac- 
tion depends quadratically on while the substrate- 
particle interaction has a linear dependence on x- There- 
fore, in the limit of small the particle-substrate in- 
teraction is the leading contribution to the total energy. 
Hence, valuable information about the model can be ex- 
tracted by simply analyzing the particle-substrate con- 
tribution to the total energy, i.e. the limiting case of 
a single colloidal particle. The potential (|3]) for a sin- 
gle particle, taken as i = 1, reads /3/7(xi, Zcoii, ^) = 
— /3/isxH(xi, Zcoii, ^)^. Figure [2j a) shows this poten- 
tial as a function of x = xi for the case of yi = and for 
a pattern of parallel stripes (as shown in Fig.ITJi), at time 
t = 0. The external field is zero and the energy has a se- 
ries of minima at the positions of the dipoles (shown in 
Fig.[2fd)), whereas the maximum of the energy is exactly 
half-way between two dipoles. Let us assume that at 
t = the particle is sitting in the energy minimum at po- 
sition xi = — 3cr, i.e. at the position of a dipole pointing 
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The particles can on average advance in discrete steps 
of 0.5 A every half period tq. We will hence quan- 
tify the transport by the (time) average current (Jx) = 
{jj" J2i ^ii^l^^iM^^ where Xi{t) and Xi{to) are the posi- 
tions of the colloidal particle i at time t and initial time 
to, respectively. With this definition, the maximum cur- 
rent measurable is (Jx) = A/tq, and an average current 
(Jx) < A/ro is an indication of a decreased efficiency of 
the transport mechanism, e.g. due to thermal motion 
or collisions among particles. Performing computer sim- 
ulations, as laid out in the next section, allows us to 
investigate the values of the external magnetic field that 
induce a particle current and the effect of the particles' 
Brownian motion as well as effects due to many particle 
interactions. 



Figure 2: (Color online) (a) Scaled energy landscape /3U as a 
function of x/a at position y — 0^ and z — ZcoW for a pattern 
of parallel stripes at time t = 0. The circle indicates the 
preferred position of the particle, (b) Same as (a) but at time 
t=To/4 in a tilted external potential. The arrow indicates 
that the particle jumps to the next energy minimum, (c) 
Same as (b) but at time t — 3to/4. (d) Magnetic field lines of 
a sequence of positive and negative point dipoles in the x — z 
plane. These dipoles generate the energy landscape shown in 
(a). 

in the positive direction. At time t = ro/4 the external 
field has positive x- and z-components. Consequently, 
it reduces the total field above the dipoles pointing in 
the positive x-direction and enhances the field above the 
dipoles pointing in the negative x-direction. Likewise 
the field between dipoles is enhanced and reduced alter- 
nately. This gives rise to an asymmetry in the energy 
landscape, as shown in Fig. [2fb). Due to the presence 
of a point of inflection in the energy, the particle moves 
towards the energy minimum that is now located above 
a negative dipole. At time t = ro/2 the external field 
vanishes, and the energy landscape is the same as shown 
in Fig. [2ja). On the other hand, at time t = 3ro/4 the 
external field has negative x- and 2;-components. The to- 
tal field above the positive dipoles is now enhanced and 
the total field above of the negative dipoles is reduced. 
The field between dipoles is again reduced and enhanced 
alternately. The point of inflection is again present in 
the energy landscape (Fig. Wi'^)) and the particle moves 
to the next energy minimum, which is now located above 
a positive dipole. It is clear from the sequence of Fig. [21 
that after one cycle tq of the external field the particle 
has covered the distance A and that the cycle can be re- 
peated indefinitely. When the inclination angle of the 
external field is such that either the z-component or the 
x-component is zero, the point of inflection in the energy 
landscape is never formed and the particle does not move 
in any preferential direction. For the point of inflection 
to form both an x- and a 2:-components of the external 
field are necessary to break the symmetry for a pattern 
of parallel stripes. 



C. Simulation Method 

The simulation box has a lateral size Lx ^ Ly — 
30 X 60 cr^ and has periodic boundary conditions in the 
X- and ^-directions. The substrate lies in the z = 
plane and contains 2600 point dipoles with dipole mo- 
ment mo = bOHoCF^ with Hq the unit of the magnetic 
field. We choose a wavelength A = 3cr, and a dipole sep- 
aration distance /S.y = 0.3a. The colloidal particles are 
constrained to move in the z = 2:coii=const. plane and 
have a susceptibility x = 0.4 arbitrarily chosen such 
that the linear term in the energy is the leading term. 
We sampled averages for 100 tq, after 5 tq of 'equilibra- 
tion' time. The long range dipole-dipole interactions are 
treated through the Ewald sum m,!!^. In order to speed 
up the interaction calculations, we pre-compute the field 
Hgub due to the substrate on a 400 x 400 grid. During 
the simulations, the field intensity is obtained by inter- 
polation of the tabulated values. The time is in units 
of the Brownian time tb = cr^ /D^ with D = ksT/^ the 
Stokes-Einstein diffusion coefficient of the particles and 
^ the friction coefficient of the solvent. In our simula- 
tions, the hydrodynamic interactions are neglected. Due 
to solvent hydrodynamics, the diffusion coefficient of par- 
ticles depends on their distance from the substrate |25l - 
32]. Therefore, simulation carried out at constant dis- 
tance Zcoii are characterized by a Brownian time tb that 
depends on ZcoW- However, this has no effects on our 
results because they are scaled by the Brownian time. 
Further many-body effects due to hydrodynamics are ne- 
glected though. 

We carry out both standard Metropolis Monte 
Carlo t33|] (MC) simulations with a small MC displace- 
ment d = O.Olcr, and Brownian Dynamics (BD) simu- 
lations [33[. The relationship between MC and BD has 
been extensively studied in the literature. Both modified 
MC scheme s llTHlQj as well as standard Metropohs MC 
simulations [2Q|-[22| give dynamical properties that can be 
in good agreement with the results of BD simulations. In 
particular, it was recently shown [22] that the dynamical 
properties obtained from Metropolis MC simulations are 
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in good agreement with those obtained from BD as long 
as the maximum step size of the MC move, is smah 
enough and the time scale in MC simulations is obtained 
according to the relations St = ad? /Qtb where a is the 
average acceptance probability of the MC moves. The 
relation was demonstrated for one particle in a arbitrary 
one-dimensional potential, and verified explicitly for a 
many-body system and various three-dimensional poten- 
tials. Scaling the time with the acceptance probability 
in equivalent to advancing the MC time only when a 
move is accepted, this concept was called 'internal clock' 
by Royall et al. [SJ. Given the presence of an oscillat- 
ing external field in our model we explicitly investigate 
the agreement between MC and BD simulations at low 
particle density. 

We carried out simulations for a set of four distances 
between the colloidal particles and the substrates, Zcow = 
1.0,1.5,1.7,2(7. For z^oW = 1.0 a the current {Jx) was 
zero in all cases while the results for ZcoW = 1-5, 1.7, 2 cr 
are discussed in detail below. 
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Figure 3: (Color online) Average current (Jx) in the direction 
perpendicular to the parallel stripes for Zcoii = 1-5, 1.7, 2.0 cr. 
(a) As a function of the period tq of the external field, (b) As 
a function of the temperature-susceptibility ratio, ksT/iUs, 
with an external field oscillation period To = IOtb- 



III. RESULTS 

First we study the low-density behavior of the sus- 
pended fluid of colloidal particles on a pattern of parallel 
lines of point dipoles, as shown in Fig. [1^. We apply a 
tilted oscillating external magnetic field, with vanishing 

component, H^^^{t) = 0, and x- and z-components 
given by H^^{t) = i^^'^ sin(27rt/ro). The current mea- 
surements were conducted using only a single colloidal 
particle. Figure [3] (a) shows the time average current 
{Jx) as a function of the oscillation period of the exter- 
nal field Hext (^) for three different values of ^^coii and with 
ksT/fis = 5 X 10~^. We find that the current is induced 
only for values of the oscillation period tq larger than a 
critical value. This result can be interpreted easily. If the 
external magnetic field is oscillating too quickly (small 
period, high frequency) on the scale of the characteris- 
tic diffusion time (Brownian time), then the particle is 
unable to follow the rapidly changing energy landscape. 
Figure [31(b) shows the x-component of the measured cur- 
rent (Jx)^ 

as a function of kBT/fig^ for three different 
values of Zcoii and with tq = IOOtb- 

We find that the external field drives a particle cur- 
rent for values of temperature over permeability, ksT/fig 
smaller than a critical value. This result suggests that if 
the thermal energy is too large, the Brownian motion 
randomizes the motion of the colloidal particles and sup- 
presses the transport. 

FigurelUa) shows the comparison between MD and BD 
simulation results for the current as a function of the pe- 
riod of oscillation tq. The two simulation techniques give 
qualitatively the same behavior including the presence 
of a critical oscillation period beyond which the trans- 
port is possible. Nevertheless, the BD predicts a critical 
period that is two times smaller than that of the MC 
simulations. We also show the period rescaled by the av- 



erage acceptance probability a as suggested by Sanz and 
Marenduzzo [22[. The procedure leads to a much better 
but not perfect comparison. On the other hand, for oscil- 
lation periods far from the critical value the comparison 
between BD and MD is very good. Figure SJ^b) shows 
that for a value tq = IOtb, both BD and MC simulations 
give the same behavior for the current as a function of 
the temperature over permeability. 

Figure [5](a) shows the value of the average current (Jx) 
as a function of the magnitude of the external field |i^ext | , 
at a fixed inclination angle Ot = arctan(i7^/i^^) = 45^ 
of the external magnetic field, ^^coii = 1.7cr and with an 
oscillation period tq = lOr^. We find that a critical 
value of the external field needs to be reached in order 
to initiate the particle transport. Figure [5fb) shows the 
value of the average current (Jx) as a function of the 
inclination angle Of at a fixed external field magnitude 
|i^ext| and for Zcoii = 1.7 a. We find that the current 
magnitude and direction can both be controlled by the 
inclination angle of the external magnetic field. Note that 
for angles Ot = ±180^, ±90^, 0^ the current is zero. These 
angles correspond to an external field with either only a 
component of the external magnetic field parallel to the 
substrate {Ot = ±180^,0^) or only a component normal 
to the substrate {Ot = ±90^). Figure [5] also shows BD 
simulation results. The comparison between MC and BD 
simulations is very good. 

We next investigate the dependence on particle den- 
sity. We characterize the system by a linear density 
p = 2aN/NxLy, where N is the total number of par- 
ticles. Here, we carried out only MC simulations. These 
posses higher computational efficiency over BD at the 
large number of particles that we are considering. Fig- 
ure [6] shows the current induced on a parallel stripes pat- 
tern, as a function of p. We find that the current remains 
roughly constant at low densities with a value slightly 
smaller than unity. This result indicates a decreased effi- 
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Figure 4: (Color online) Comparison between MC, MC scaled 
with the acceptance probability a, and BD simulation results 
for the measured average current {Jx) for Zcoii = 1.5 a. (a) 
As a function of the period tq of the external field, (b) As 
a function of the temperature-susceptibility ratio, ksT/jjis^ 
with an external field oscillation period tq = lOrs- 



Figure 6: (Color online) Average current {Jx) for a pattern of 
parallel stripes, Zcoii = 1.7 a, ro = lOrs and ksT / iis — 0.05 
as a function of the linear density p. 
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Figure 5: (Color online) Average current {Jx) in the direction 
perpendicular to the parallel stripes for Zcoii = 1.7 a. (a) As 
a function of the magnitude |ifext| of the external field at 
an inclination angle Ot — aYctan{H^ /H^) = 45^. (b) As a 
function of the inclination angle Ot for a magnitude |ifext| = 
50 mo/cr^. 



ciency of the transport mechanism due to the presence of 
other particles. Interestingly, the current decays to zero 
at densities larger than p :^ 2. At this density, all stripes 
are filled with particles. If more particles are present, 
they must fill positions that are not ideal for the trans- 
port mechanism. These excess particles effectively jam 
the transport. 

We next analyze the motion of the particles suspended 
over a zigzag pattern with an external field perpendic- 
ular to the substrate. The x- and ^/-components of the 
external field vanish, H^^{t) = 0, and the perpendicu- 
lar z-component, H^^^{t) = sin(27rt/ro), is oscillating 
in time with period tq. Analysis of the particles tra- 
jectories shows that transport is achievable only in the 
vertex region of the zigzag. Therefore, to reliably mea- 



sure the current, we initialized the simulations with a 
single colloid randomly positioned along the x direction 
but localized at the vertex. Figure [Tfa) shows the cur- 
rent {Jx) as a function of the period of oscillation tq, for 
three different values of Zcoii and for 6^ = 68^. As for 
the case of parallel lines, the period of oscillation of the 
external field needs to be large enough in order to induce 
a current. Furthermore, we find that for Zcoii = 1.5, 1.7<j 
the current is negative, while for Zcoii = 2.0cr the current 
is positive. Figure [7Jb) shows the average current as a 
function of the zigzag angle 0^ for external field Hz = 10 
mo/ (7^ and tq = lOOr^. The current direction is differ- 
ent for different distances Zcoii- In particular it is positive 
for Zcoii = 2. Oct, and negative for ^^coii = 1.5cr. Interest- 
ingly, for the intermediate value Zcoii = 1.7cr, the current 
is negative for large angles, but positive at smaller zigzag 
angles. The current vanishes for Oz = 90^, corresponding 
to the case of parallel stripes. We stress that the current 
is only obtained at the zigzag vertex because the exter- 
nal field only has a normal component. This means that 
the local relative position of the dipole located around 
the vertex plays a fundamental role in the creation of 
the point of inflection in the energy, that is necessary for 
the transport mechanism (see Sec. IIIBj) . Nevertheless, if 
a tilted external field is applied to the zigzag patterns, 
transport is achieved also away form the vertex. Fig- 
ure [8ja) and (b) show the comparison between MC and 
BD simulations for Zcoii = 1.5cr. We find a discrepancy 
in the predicted value of the critical period. Contrary to 
the case of the parallel lines pattern, the BD predicts a 
critical period that is two time larger than MD. Rescaling 
with the average acceptance probability leads to a worse 
comparison. On the other hand, the comparison is very 
good for the current as a function of the zigzag angle for 
a value of the oscillation period tq = Ttb far from the 
critical value. 
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Figure 7: (Color online) Average current {Jx) for a zigzag 
pattern with ksT / fis = 5 x 10""^ for Zcow — 1-5, 1.7, 2.0a. (a) 
As a function of the oscillation period tq of the external field 
for Oz = 68^. (b) As a function of the zigzag angle Oz for 
To = 7tb. 
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Figure 8: (Color online) Comparison between MC, MC scaled 
with the acceptance probability a, and BD simulation for 
the measured average current (Jx) for Zcou — 1.5 a. a) As 
a function of the oscillation period tq of the external field 
for Oz = 68^. (b) As a function of the zigzag angle Oz for 
To - 7tb. 



IV. CONCLUSIONS 

We have studied a simple model for the transport of 
colloidal particles suspended at a fixed distance over a 
magnetic patterned substrate with MC and BD computer 
simulations Magnetic dipoles were distributed in two spe- 
cific patterns, namely parallel stripes and zigzag stripes. 
We analyzed the effect of an oscillating external magnetic 
field that was applied to the system. 

For the case of parallel stripes we found that the cur- 
rent magnitude and direction was controlled by the tilt 
angle of the external field and that the effect was re- 
liably obtained in a wide range of ratios between tem- 
perature and solvent permeability. Furthermore, a net 
current was measured only when the period of oscilla- 
tion was greater than a critical value. For the case of 



zigzag stripes a current was obtained using an oscillat- 
ing external field normal to the substrate. In this case, 
transport was only possible in a small region of the pat- 
terned substrate, namely near the vertex of the zigzag. 
This result opens up the possibility to transport colloidal 
particles in a very narrow stream. Furthermore, the cur- 
rent magnitude and direction was found to be controlled 
by a combination of the zigzag angle and the distance of 
the colloids from the substrate. The comparison between 
MD and BD is overall qualitatively very good. We find 
quantitative agreement for values of the period of oscil- 
lation of the external field far from the critical period, 
while the two simulations techniques predicts quantita- 
tively different values of the critical period beyond which 
transport of particles is possible. 

The mechanism behind the transport of the colloidal 
particles is a consequence of the changing energy land- 
scape. The sum of the oscillating external field and of the 
substrate's magnetic field results in an energy landscape 
that changes in time. The Brownian motion enables the 
particles to locally sample the phase space and follow the 
energy landscape towards the local (in space and time) 
energy minimum. Colloidal transport is hence achieved 
when the particles are able to "follow" this landscape. 
The mechanism explained here is the same as the one 
described by Dhar et al. [l^ as a deterministic ratchet. 
Yellen et al. [l2| found the same mechanism and describe 
it as particles following a traveling wave. That is, the par- 
ticle is transported by the translating infiection point in 
the energy landscape. We find that transport of the para- 
magnetic colloidal particles is possible for a large set of 
model parameters. The magnetic patterns can be created 
by deposition of discrete magnetic islands [lo|, fl 2| . [l 3|. 
Despite that the field generated by a garnet film |l5j is 
quantitatively different from the one produced by an ar- 
ray of discrete dipoles, the differences are surprisingly 
small (see appendix [A|) . Therefore, we expect that the 
behavior for the transport of particles on top of garnet 
films is similar to the one shown by our model. 

Controlling the deposition pattern means controlling 
the behavior of the nano or micro magnetic particles. 
As a consequence lab-on-chip devices with well defined 
functions can be envisioned. Our model and method 
can be easily applied to different and more complicated 
patterns, like for example a combination of parallel and 
zigzag stripes. Other possible extensions of the current 
work include the study of substrate boundary effects. In 
this work we applied periodic boundary condition, but it 
would be interesting to study the transport of colloidal 
particles on top of finite discrete patterns. 

Furthermore, given the range of phenomena shown by 
two-dimensional colloidal suspensions of paramagnetic 
particles trapped at a liquid-air interface (see for exam- 
ple Ebert et al. [35| and references therein) it would be 
interesting to explore in more detail the effect of a pat- 
terned magnetic substrate on the phase behavior and the 
dynamical properties of two-dimensional fiuids. 
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Figure 9: (Color online) (a) Comparison between the energy 
of one paramagnetic particle in the magnetic field of a garnet 
film and the dipole array, (b) Like (a) but zoomed in to 
highlight the differences. 
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Appendix A: Garnet Film 

Tierno et al. [15] calculated the magnetic field above 
the garnet film, for a pattern of parallel stripes aligned 
in the y direction as H = VRe[<l>] with the potential 

^ = — dilog(l — exp(-— + Im(/i))) 

z 4z7r 
- -dilog(l + ex.p{—w - lm{h))), (Al) 

TT 2 

with w = X -\-iz and h = Hext • {^x — '^^z) and the dilog- 
arithm function dilog(t) = Ji dt-^^- Figure [9] shows 
the comparison between the energy of the garnet film 
from the potential ()Aip and the energy of an array of 
dipoles from equation (|3]). The energies were shifted by 
the mean value (E) and rescaled by the maximum value 
of the energy £^max- Both with and without external field 
the differences between the two energies are small, val- 
idating the use of the set of discrete dipoles as a good 
approximation for patterns on garnet films. 
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